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1 Introduction 

Propagation and electron cyclotron resonance (ECR) damping of electron Bernstein 
waves (EBWs) in spherical tokamaks (ST) is usually analyzed assuming that absolute 
value of the tokamak magnetic field \B\ increases inward the plasma. In this case the 
perpendicular index of refraction n± grows as the wave approaches the ECR layer, 
where the wave is fully absorbed regardless of the resonance harmonic number PJ. 
However, in addition to traditional regimes with monotonously increasing \B\, regimes 




Figure 1: Total magnetic field B profile in NSTX tokamak mid-plane 0. 

with "magnetic wells" also occur in STs [2]. The the magnetic field profile inversion 
modifies significantly the whole picture of the wave propagation and damping. Since 
the magnetic wells may become quite common with further improvement of ST per- 
formance (fig.l), analysis of such configurations is of interest for assessment of EBW 
plasma heating an CD perspectives. 

Consider CMA diagram (fig. 2) illustrating why this case was not considered in 
literature yet. In the conventional tokamaks EBWs produced via the linear con- 
version of the incident electromagnetic waves close to UHR are capable for prop- 
agating to the fundamental harmonic along dotted blue arrow. If the inequality 
d[ln{ne)]/dx > 2u'^^/Up^d[ln{\B\)]/ dx is fulfilled (it's the case of STs) EBWs can propa- 
gate along solid red arrow to the nearest EC harmonic. In this paper we consider basic 
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Figure 2: CM A diagram. In the case of the conventional tokamak EBWs can propa- 
gate from UHR to the fundamental harmonic (dotted arrow). In the ST EBWs can 
propagate to the second EC harmonic (solid arrow). This case will be considered in 
this paper. 



features of the EBWs propagation and damping for the second cyclotron harmonic, 
which is now the lowest possible resonance harmonic, in a slab model. This case is 
illustrated in fig.l by solid red arrow. 



2 Bernstein waves in electrostatic approximation 

Assume all plasma parameters depending on the single dimensionless co-ordinate x 
scaled in units of c/uj with uj being the wave frequency. Suppose that the magnetic 
field B is along the z - axis and both v = oj^J'jJ^ and q = uj/uce grow inward the 
plasma. 

We begin analysis of the EBW behavior assuming validity of the electrostatic approx- 
imation. Then the dispersion relation is 

e (n_L,n||,x) = 0, (1) 

where e = {niSikUk) jr? is the longitudinal dielectric function, - the plasma dielectric 
tensor elements, n_L = kxc/uj and n\\ = k\\c/uj. Since the UHR is known to be the cutoff 
for electrostatic EBWs, one expand the dielectric tensor components eij,i,j = {x,y,z) 
in power of A = k\pl = {nj_qf3)'^ /2, where pe is the electron Larmor radius. As it will 
be seen, n± remains within the limits 1 <^ n_L < f3~^, where f3 = I'te/c, ^te = y^2Te/me, 
in the whole region between the UHR and the q = 2 resonance, so that A < 1 here. 
For qualitative investigation we keep only zero- and first-order terms. Then 



where e^'^^ — (s^^I'iA. + ^iz^f) /"-^ 

.2 

g2 - 1 



4^2 = i--?^,4^i = i-- (3) 



arc the dielectric tensor elements in the cold plasma. The electron Bernstein waves are 
produced via the linear conversion of incident electromagnetic waves with ny < 1, while 
applicability of the electrostatic approximation requires that n± ^ 1. This permits 
one to omit terms proportional to n| and n| in Eq.(2) reducing the dispersion relation 
(l)to 

^"1^0. (4) 



The solution to this equation 



where s^^^ < and 1 < q < 2 shows that in the electrostatic approximation the EBW is 
confined between two cut-offs. One of them is the UHR, for which €l^].{x) = 0, and the 
other one is the cyclotron resonance q = 2. The validity condition for the electrostatic 
approximation is ^ \sik\- We consider here typical for STs high - density plasmas 
with V ^ 1, then the validity condition becomes 

> V. (6) 

This condition breaks down close to the UHR and ECR. We do not consider first 
of these regions because analysis of wave behavior there is the subject of the mode 
coupling theory In the dense plasma outside immediate UHR vicinity, Eq. (5) becomes 



so that the characteristic value of n± is f3 ^. In the 2"^* ECR vicinity, where the 
inequality (6) is broken down a full - wave treatment is required. 



3 Full-wave equation near q = 2 resonance 

To obtain a traceable full - wave hot plasma dispersion relation for EBWs with n± ~ 
and n_L ^ near the q = 2 resonance we make some simplifications. Consider first 
the dielectric tensor elements Sik- It is well known that the elements can be presented 
as an infinite sums over cyclotron harmonic number s, — oo <^ s <^ oo with each term 
of the sum related to the ECR q — s. Since parameter A is small, we calculate the 
resonance (s = 2) terms of e^k up to the first order in A using the zero order (cold 
plasma) approximation for non - resonant terms. In this approximation, the elements 

^xzi ^zxi ^yzi ^zy Vanish, Sxx — ^yy — ^-L? ^xy — ^yx — and 

2n\\p 2n||/5 
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(8) 



here 

= -v\q\/ {q^ - 1), Z is the plasma dispersion function defined according to 
jSj with the argument C, = {q — 2)/qn\\(3. Consider now the wave equations 



[n^]^ — £:_L^ Ex — igEy — n\\n±Ez = 

tgEx + [n^ - 5x) Ey = (9) 
-n\\n^Ex + (nj - e^z) E^ = 

Absolute values of dielectric tensor elements here are of order v ^ 1 and, therefore 
Ez <^ Ex- Omitting Ez in the first equation, that means neglecting terms of order 
n±/v and 1/v compared to unity, obtain the dispersion relation 

nle± - {ei_ - g) {e± + g) =0. 

Using Eq.(8), we obtain at t; ^ 1 

„.,,_4¥^^_!4qi^M.o. (10, 

Close to UHR, at large nj_, n\ ^ v, the electrostatic approach is valid, the first term 
in (10) dominates and we returns to the electrostatic equation (5). In the ECR vicinity, 
at small n±, n± ^ v, this term can be omitted and the solution is found explicitly: 



^ /?^z((g-2)/gn||/3)- 

Thus, in the vicinity of the q = 2 resonance, where n±{x) goes down with x, there is a 
spatial region v"^ > n']_> v where Eqs.(5) and (11) are valid simultaneously. 
The dispersion relation (10) can be also used in the "relativistic" case riy < (3 if the 
"non-relativistic" plasma dispersion function Z is replaced by a proper relativistic 
dispersion function. 



4 Wave behavior in the 2 ECR layer 

Close to the 2"°' ECR layer we consider, for simplicity, as x dependent only the resonant 
factor (g — 2)/g in the argument of the Z function, putting {q — 2)/q = x/l, where / is 
the dimensionless characteristic scale-length of the magnetic field variation. Using the 
asymptotic expression Z ~ —1/^, where ^ = x/L and L = l(3n\\, at |^| ^ 1, one find 
that at negative x outside the resonance layer Eq.(ll) differs from the electrostatic 
equation (5) only by a constant factor in the right-hand side. Inside the EC layer 
at 1^1 ~ 1 the function Z is a complex function with \Z\ ~ 1 and Re{Z) ~ Im{Z). 

Solution to Eq.(ll) is also complex with \n±\ ~ (n\\/f3^ , so that the applicability 
condition for the WKB theory in the resonance region takes the form 

n^L = > 1, (12) 
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The estimation for the relativistic case can be obtained from this equation by putting 
n\\ ~ (3. We assume that the WKB approximation is apphcable in the region outside 
the resonance layer where |nj_| grows. Suppose now that inequahty (12) is satisfied and 
consider solution to Eq.(ll) inside the resonance layer. This solution in a parameter 
free form 

N = ±Z-i/^ N = n^ ((3/{An{)f^^ (13) 
is shown in fig. 3 on the plane {Re{N), Im{N)). Brunches located in the upper 
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Figure 3: Plane {Re{N) , Im{N)) . Left branch corresponds to incident wave, right one 
- to reflected. In WKB approximation branches are not intersected. Reflection is the 
result of WKB approximation breaking. 

and lower complex A^-plane correspond to the signs (— ) and (+) in Eq.(13), re- 
spectively. Consider first the brunch = — The fragment of the curve 
with Re{N) < 0, Re{N) ^ Im(N) belongs to the region x < 0, \x\ ^ L 

,1/2 



- (4|x|n||/(/?/))' 



At 



when the asymptotic expression for Z can be used and n± 
Since {Im{Z)) < 0, the whole curve lies in the upper complex plane 

1 /2 

|x| ^ L we have n± ^ i (^4|x|n||/(/5/)^ . The second curve describes transition from 

n± (^4|x|?T,||/(/5Z)^ ^ to n_|_ — > —i (^4:\x\n\\/ {(31)'^ ^ via the lower complex A^ plane. 
Note that the waves under consideration have their phase and group velocities direct- 
ed oppositely. Therefore, the incident wave has asymptotically a negative n± in the 
propagation region x < 0. This wave is described by the A^ = —Z^^^"^ dispersion curve. 
At large positive x this mode vanishes exponentially. The second brunch describes the 
wave reflected from the ECR. This mode grows exponentially in the non propagation 
region |a;| 3> L. Three important conclusions can be drawn from this analysis. First, 
incoming waves incident on the ECR layer from the high-field side are not converted 
in the resonance region into outgoing EBWs with large n± propagating on the low 
field side of the ECR. Instead, the incident waves become non-propagating beyond the 
resonance layer. Second, two effects are simultaneously responsible for Im{n±): the 



ECR damping and the wave transition into the non-propagation region. These two 
contribution can not be separated. Finally, since the function has no zeros 

at finite |^|, two brunches of the dispersion curves are separated in the whole complex 
X plane. As a result, refiection from the ECR layer can only be due to approximate 
nature of the WKB theory. 



5 Reflection from ECR layer 

We analyze reflection of EBWs from the q — 2 cyclotron resonance with the use of 
the model wave equation obtained from the dispersion relation (13) by replacement 
n± — > —id/dx: 

U"{x) +n±{xfU{x) = 0. (14) 

Here is given by the right - hand side of Eq.(ll) and wc assume that function 
U is related linearly to the wave field components. More definite interpretation of 
this function is not required. In the WKB approximation Eq.(14) has two hnearly 
independent solutions 

Ui = ^y^exp ^—i J n±dx^, U2 = -^^exp J n±dx^ (15) 



where a is an arbitrary constant. As it has been mentioned before, Ui and U2 represent 
incoming and outgoing waves, respectively. The solution to Eq.(14) describing the ECR 
damping and reflection of the wave incident from the high- field side vanishes at |x| ^ L 
and has the asymptotic form 

U = U^ + RU2 (16) 

at a; < 0,|,t| ^ L. Here i? is a constant and the quantity \R\'^ is the energy reflection 
coefficient. Presenting Eq.(14) in the form 

U"{x) + (Vo + Vi)U(x)^0, (17) 

where Vq = n\ — Vi^ Vi = {n^^^'^\ nj^. One can easily see that functions are 



exact solutions to the equation + ^of^i,2 = and V\ is small, of order (|n^L|)~^ 
compared to Vq, if the WKB method appficability condition \n\\ L 1 is satisfied. In 
this case the reflection coefficient R can be calculating by treating the term in Eq.(17) 
as a small perturbation and using C/i as a zero order solution: 



or explicitly. 



R=\\ VxUldx (18) 



where the prime denotes differentiation with respect to and Fq = 2//3^/^ny'^^. For 
Eq.(19) validity, Fq must be large compared to unity. Since the integrand has no 
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saddle points at finite |^|, and contribution to the integral is spread over the whole 
|.^| < 1, the possibility of obtaining an accurate analytical estimation of the integral 
seems rather doubtful. Numerical evaluation shows that \R\ is negligibly small at 

ro>i. 

In the limiting case \k±\L < 1 opposite to the WKB one, the reflection coefficient can 
be found by presenting Z^^(^) ~ — ^ + V and treating V as a. perturbation. Substi- 
tuting these expressions into Eq.(14) and introducing the new independent variable 

T = kox {P'^l/iy^^^, obtain 

U" - (r - 7^ (r/7)) U = 0, (20) 

where 7 = Tq^^ is a small parameter. In the present case the reflection coefficient is 
close to unity. Calculating the correction to it in the lowest order in 7 yields 

1^1 (21) 
where = iTAi{0) = 1.12. The comparison of the reflection coefficient, R, calculated 
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Figure 4: Absolute value of the reflection coefficient \R\ for the non-relativistic (ny ^ j3) 
case. 

numerically and its analytical estimation (21) is presented in fig. 4. The reasonable 
agreement is demonstrated. Note in conclusion that standard ray tracing method 
using the full wave dispersion relation or the approximate dispersion relation (11) can 
be used for high-field side propagation only if the wave is damped in the periphery of 
the ECR, where Im{n±) <C Re{n±). Otherwise the wave penetrates into the region of 
strong cyclotron damping where the ray tracing method is inapplicable. 
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6 Summary 



Existence of EBWs in the region between the UHR and ECR requires inhomogeneous 
plasma density. 

Waves in this region are adequately described by the approximate full-wave dispersion 
relation (10). Incoming waves incident on the ECR layer from the high-field side are 

not converted in the resonance region into outgoing EBWs propagating on the low 
field side of the ECR. Instead, the incident waves become non-propagating beyond the 
resonance layer. 

Decreasing of the wave amplitude within ECR layer is due to combine effect of the 
ECR damping and non-propagation. 

In the WKB approximation, the waves are fully damped in the ECR layer. Reflection 
from the ECR layer is only due to approximate nature of the WKB theory. 
Standard ray tracing method can be used for high-field side propagation only if the 
wave is damped in the periphery of the ECR layer. Otherwise the wave penetrates in- 
to the region of strong cyclotron damping where the ray tracing method is inapplicable. 

Acknowledgment: This work has been supported by RFBR 04-02-16404, 02-02- 
17683, Scientific School grant 2159.2003.2. 



References 

[1] A D PiLiYA, A Yu Popov and E N Tregubova, Plasma Phys. Control. 
Fusion 45 (2003) 1309 

[2] B. Jones, G. Taylor, P. C. Efthimion, and T. Munsat, Physics of Plasmas 
11 (2004) 1028 

[3] D.G. SWANSON, Mod. Rev. Phys. 67 (1995) 837 



8 



